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1 Introduction 



The scalar field model with potential U(4>) = ^\<\> 2 + ^4> 4 + f,0 6 is the simplest model 
exhibiting a rich phase structure and for studying tricritical phenomena in both two and 
three dimensional systems [see, e.g. Refs. and ||]. In terms of the phase diagram 
for the model, a tricritical point can emerge whenever we have three phases coexisting 
simultaneously. For the above potential, in the absence of corrections due to fluctuations, 
one can have a second order transition in m when the scalar field mass vanishes and 
A > 0, a > 0. A first order transition happens for the case of A < 0, a > 0. The tricritical 
point occurs when the quartic coupling constant vanishes (with m = 0, cr > 0). A study 
of the phase diagram for this model in 3D, at finite temperature, was recently done 
and it has been shown that there is a temperature /3~ l (m, A, cr) for which the physical 
thermal mass mp and coupling constant \p vanishes, thus characterizing the tricritical 
point. 

There are many possible applications associated with the model with potential U((f>) 
with 6 interaction. For instance, in D = 2, it is known H that the minimal conformal 
quantum field theory, with central charge 7/10 (the tricritical Ising model) is in the 
same universality class, in the scaling region near the tricritical point, of the Landau- 
Ginzburg model with the above potential. This (ft 6 potential model can then be thought 
of as the continuum realization of the Ising model with possible applications in, e.g., 
the description of adsorbed helium on krypton-plated graphite ||, in understanding the 
statistical mechanics of binary mixtures, such as He 3 — He 4 etc. These are just a 
few examples of systems exhibiting tricritical phenomena in condensed matter physics. 

In field theory in general, the 6 model has been used in the study of polarons, or 
solitonic like field configurations on systems of low dimensionality 0. Also, a gauged 
(SU(2)) version of the 6 potential model in Euclidean 3D, has recently been used to study 
a possible existence of a tricritical point in Higgs models at high temperatures |7[]. In this 
context the tricritical point is characterized by the ratio of the quartic coupling constant 
and the gauge coupling constant of the effective three-dimensional theory, obtained from 
the 3 + 1 D high temperature, dimensionally reduced SU(2) Higgs model. This can be 
particularly useful in the context of the study of the electroweak phase transition. 
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There are then many reasons that make the 6 potential model an interesting model 
to be studied. In this paper, we will be particularly interested in studying the regime of 
parameters for which: 



m > , A<0, a > and 



>0. (1) 



In this case U{<fi) has three relative minima <pt± and 0/ (see Fig. 1). For these parameters 
the system has metastable vacuum states and it may exhibit a first order phase transition. 
The states of the classical field theory for which = <p t ± are the unique classical states 
of lowest energy (true vacuum) and, at least in perturbation theory, they correspond to 
the unique vacuum states of the quantum theory. The state of the classical field theory 
for which = 0j is a stable classical equilibrium state. However, it is rendered unstable 
by quantum effects, i.e., barrier penetration (or over the barrier thermal fluctuations, at 
finite temperatures). 0/ is the false vacuum (the metastable state). 

We will compute the vacuum decay (tunneling) rate at both zero and finite temper- 
ature. For calculation reasons we will restrict ourselves to the thin wall approximation 
for the true vacuum bubble (or bounce solution). We then consider the energy- density 
difference between the true and false vacuum as very small as compared with the height 
of the barrier of the U(<f)) potential |], |TB[. From this, we are able to give the explicit 
expression for the bounce and also to qualitatively describe the eigensolutions for the 
bounce configuration. The paper is organized as follows: In Sec. II we obtain the bounce 
field configuration for the model and we compute the Euclidean action in the thin wall 
approximation. In Sec. Ill we calculate the vacuum decay rate at zero temperature and 
we discuss the eigenvalue equations obtained for the bounce configuration within the ap- 
proximations we have taken for the bounce. In Sec. IV we compute the nucleation rate 



at finite temperature, following the procedure given in [14]. In Sec. V the concluding 
remarks are given. In this paper we use H = c = = 1. 
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2 The Vacuum Decay Rate and the Bounce Solution 

Let us consider a scalar field model, in three-dimensional space-time, with Euclidean 
action given by 



with 



S E W = Jd s x E [i (8^f + (7(0)] , (2) 



As discussed in the introduction, we are interested in the regime for which the parameters 
in (|3]) satisfies (P, such that the potential U((j)) exhibits non-degenerate local minima and 
then metastable states. The picture we have in mind is that once the system is prepared 
in the false vacuum state, it will evolve to the true vacuum state by tunneling (at zero 
temperature) or by bubble nucleation, triggered by thermal fluctuations over the potential 
barrier. 

In the case of quantum field theory at zero temperature the study of the decay of 
false vacuum was initiated by Voloshin, Kobsarev and Okun || and later by Callan and 
Coleman H |10| , who developed the so called bounce method for the theory of quantum 



decay. In this context the decay rate per unit space-time volume V3 is given by 



T_ _ (AS E \ 3/2 

v 3 ~ v^r 



det'i-VE + U"^)) 
det (-D E + U»{(f> f )) 



-1/2 



e- AS *(l + 0(h)) , (4) 



where U E = £ + + ^, U"{<f>) = ^ and AS E = S E (<p b ) - S E {<f> f ). S E (<p b ) is 
the Euclidean action evaluated at its extreme (specifically a saddle point), <fi = fa, where 
4>b is the bounce: a solution of the field equation of motion, SS E /S(j)\ ( f )=( f >b = 0, with the 
appropriate boundary conditions. The prime in the determinantal prefactor in (^) means 
that the three zero eigenvalues (the translational modes) of the [— O e + U"(<j)b)} operator 
has been removed from it. 

As it was shown by Coleman, Glaser and Martin |]TTJ, the solution that minimizes S E 
is a spherical symmetric solution, r 2 = r 2 + x 2 + y 2 (in D dimensions the solution has 
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0(D) symmetry) and then can be written as the solution of the radial equation of 
motion 



d 



, b 2 d<f) b 
dr 2 r dr 



(5) 



with the boundary conditions: lim^oo 0&(r) = <pf and -P- 



dr l r=0 



At finite temperature the calculation of the decay rate were first considered in ]T2| 



in the context of quantum mechanics and later by Linde []13||, for quantum field theory. 



In Ul3[ , it is argued that temperature corrections to the nucleation rate are obtained 
recalling that finite temperature field theory (at sufficiently high temperatures) in D = 
d + 1 dimensions is equivalent to d-dimensional Euclidean quantum field theory with h 
substituted by T. At finite temperature, the bounce, 4>b = 4>b{p) (p = |r|), is a static 
solution of the field equation of motion: 



d\ 



(6) 



dp 2 p dp 

with boundary conditions: lim^oo 4>b(p) = 0/ and ^-\ p=0 = 0. The vacuum decay 
rate, or bubble nucleation rate in this case, is proportional to e~ AE ^ T , where AE is the 
nucleation barrier, given by (f3 = 1/T) 



AE 
~T~ 



(3 J rf 2 x [C E ((j) B )-C E {(j)f)\ 



(7) 



where £g is the Euclidean Lagrangian density. 

The problem of the computation of the nucleation rate at finite temperature was 
recently reconsidered by Gleiser, Marques and Ramos , who have used the early works 



of Langer [15]. In this context the nucleation rate is given by 



LELL |"det [-D E + U" {4>b% 

lm — 

7T [det[-D B + ?7"(0 / )] / 3 



-1/2 



e~ AE / T (1 + 0{h)) , 



(8) 



where E 2 _ is the negative eigenvalue associated to the operator [—D E + U"(<Pb))- 

To compute the vacuum decay rates, we need, therefore, to solve Eq. (|5|) for 0&(r) 
(at zero temperature), or Eq. (|6]) for 4>b(p) (at finite temperature). In the model studied 
here, we can give an approximate analytical treatment for the bounce solution in the so 
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called thin-wall approximation, in which case the energy-density difference between the 
false and true vacuum states can be considered very small, as compared to the height of 
the potential barrier: e = U(<fif) — U(<fi t ) <C U{4>2) (see Fig. 1). By interpreting (fib as 
a position and r as the time, Eq. (|5[) can be seen as a classical equation of motion for a 
particle moving in a potential —U((fi) and subject to a viscous like damping force, with 
stokes's law coefficient inversely proportional to the time. The particle motion can then 
be interpreted as it was released at rest, at time zero (because of the boundary condition 
^\ r =o = 0). From (§) it follows that 



d 
dr 



dr 



u(<p b 



dr 



< 



(9) 



meaning that the particle loses energy. Then, if the initial position of the particle is chosen 
to be at the left of some value <fi = (fii (see Fig. 2), it will never reach the position (fif. 
Now, for (fib{r) very close to (fit, we can linearize Eq. ([5]) to obtain 



dr 2 



2 d 
r dr 



(10) 



where fi 2 = U((fi t ). The solution of (|T0| ) is 



(0(0) 



^sinh(/xr) 
fir 



11) 



Therefore, if we choose the position of the particle to be initially sufficiently close to (fit, we 
can arrange for it to stay arbitrarily close to (fit for arbitrarily large r. But for sufficiently 
large r, r = R, the viscous damping force can be neglected, since its coefficient is inversely 
proportional to r. And if the viscous damping is neglected, the particle will reach the 
position <fif at a finite time R + AR. Then, by continuity, there must be an initial position 
between <fi t and (fix for which the particle will come to rest at <fif, after an infinity time. 
From the above arguments, we can find a general expression for (fib, valid for e <C U(4>2)- 
In order to not lose too much energy, we must choose 0b(O), the initial position of the 
particle, very close to <fi t . The particle then stays close to (fit until some very large time, 
r = R. Near R (between R — AR and R + AR, AR <ti R) the particle moves quickly 
(according to Eq. (|5|), neglecting the viscous damping force), through the valley in Fig. 
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2 and it slowly comes to rest at <f)f, after an infinity time. Thus, we can write for the 
bounce the following expression ((f) f = 0) 

iH < r < R - AR 

Mr) = { <j> wa n{r-R) R-AR<r<R + AR , (12) 

fi f R + AR < r < oo 
where wa n satisfies the equation 



d 



2, 



Vail 



t/'(0wall) • (13) 

ar z 

Eq. (0) is the thin- wall approximation for the bounce solution. From Eq. (0), we 
obtain 



o 



— — = r . 14 

217(0 



By rewriting U (4>) as 



6! 

with 



U{<l>) = T MP-ti) -772, ( 15 ) 



, 2 16! A 



and 



7= *i 

r 4 



6!/AV 2 
- — - 2m 2 
a I 4! ) 



;i7) 



then, by neglecting in (|15[) the term proportional to 7 (valid in the thin-wall approxima- 
tion), we obtain that 



/* 

Jo 



'2§<f> (0 2 - 



r . 



The above integral is straightforward and the solution for wa ii(r) can be written as 



7 



0wall( r ) 



1 + exp (^4>li 
This solution is shown in Fig. 3. 

Using (|19D , we obtain for ASe the expression 



(19) 



AS, 



47r / drr 2 
Jo 



5 * + W ' 



rR-AR poo 

4tt / drr 2 U((j) t ) + 4vr / drr 2 £/(0) + 
Jo Jr+ar 

^wall \ i rr/j. \ 



i-R+AR 

47r / (irr 2 
Jr-ar 



dr 



Since Ai? <C i?, in the last integral of (|20"D we can take r k, R and we obtain 



(20) 



AS E ^ ~^e R 3 + 47tR 2 S , 
with <So, the bounce surface energy density, given by 



(21) 



d(/>yj2U(<t>) « /° #7^0 (0 2 - 0o) 



r<f>f , rO 

dfhJ2TI(fh'\ sa / 

/0 O V 6! 

where we have neglected in U(<f>) the term 7</> 2 /0q and we also used 
the above integral, we obtain for Sq the result: 



(22) 
Evaluating 



5 » = tV 2 6! 



(23) 



In the next two sections, we deal with the evaluation of the determinantal prefactor 
appearing in Eqs. (£|) and (|8|) and we obtain the subsequent radiative (1-loop) corrections 
to (1211). 



3 The Vacuum Decay Rate at T = 

Let us consider, initially, the eigenvalue equations for the differential operators appearing 
in (©: 
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and 



[-D E + U"{<t> b )\ il> h {i) = E 2 b (i)Mi) 



(24) 



[-U E + U"{<P s )]^f{j) = E%j)^ f {j) . 
We then have for the determinantal prefactor of (^), the following 



(25) 



K 



det' [-U E + U" ((/>„)] 



-1/2 



det [-D E + U"((j) f )} 

"det' [-n E + u"(<fo)) 



exp 



exp 



exp 



2 



In 



det [-n E + U"((f) f )} 

'2 
b 



n< % 2 (* 



L n^|(j) 
E^i^ 2 «l-E ln ^ 2 (i) 



(26) 



Since the bounce can be approximated by a constant field configuration for r < (R 
AR) fa R, we can write for K, in the thin wall approximation, the expression: 



K 



-TTi? 3 

3 



d 3 p 
(27p 



In 



2 f 

walU 



5>^(j) 



(27) 



where E 2 ^{p) = p 2 + U" (<f>t(f)) • The integral in (|27|) can then be identified as the one 
loop correction to the classical potential, while the remaining terms represent the quantum 
corrections due to fluctuations around the bounce wall |l6fl . Then, by using Eqs. ( pT|) 
and (p7|) in Eq. (f|) and following |9|, pTOf] , we obtain 



r 

v 3 



V 2tt 



3/2 



exp 



; nR 3 AU e f f - 4nR 2 (S + <S a 



(2? 



where <Si, is the term giving the 1-loop quantum corrections to fluctuations around the 
bounce wall, 
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Si 



AnR 2 



(29) 



where £7 2 all (z) are the eigenvalues of —D E + £/"(0 wa ii(r — R). In ([2^) we have also that 
AU e fi = U e e((f>t) — U c ff((j)f), where U e g is the one loop effective potential, given by [[L7| 



M0) = W) + ~/|^ln 



p 2 + U"{<t>) 



(30) 



p2 _|_ m 2 

The ultraviolet divergence in ([3(J) can be handled in the usual way. Integrating over p Q 
and by using an ultraviolet cut-off, A, for the space-momentum, we obtain 



u cS (<p) = u{4>) - - 



2 , 2^/ 2 
p + m 



3/2 



Using that (1 + xf' 2 = 1 + 3/2x + 3/8x 2 



., we obtain for U e s the expression 



(31) 



Uestf) = U{(f>) 



1 



12 

+0(1/A) 



™ 2 + \<P 2 + ^ 



3/2 



- m 3 - |y(A) + 



(32) 



where Y(A) = A 2 + m 2 . The divergent terms in ([52) are proportional to 2 and 4> 4 but 
not to 6 . Then, only the mass m and A need to be renormalized (a r = a). From the 
usual definition of renormalized mass m r and coupling constant A r , 



d 2 U eS {4>) 



if>=0 



and 



A, 



d*UeB(<l>) 



0=0 



(33) 



(34) 



and writing the unrenormalized parameters in terms of renormalized ones, we obtain for 
the renormalized one loop effective potential, the following expression: 



M0) = UM) - ^ 



3/2 



rrr r - -m r A r 0^ 



3A 



r i 4 



32m r 



16 Y 



(35) 
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where U r means the tree level potential expressed in terms of the renormalized quantities. 
For convenience, from now on we drop the r subscript from the expressions and it is to 
be understood that the parameters m, A and a are the renormalized ones, instead of the 
bare ones. Then, AU e g can be written as (since 4>f = 0) 



off 



3ma 
"IT 



(36) 



We now turn to the problem of evaluating the eigenvalues E^^ii) of — O E +U" [4> wai n(r — 



R)], which appears in (27). This is not an easy task. In fact, only in a very few examples 
this has known analytical solutions, as, for example, for the kink solution in the (A0 4 )d=2 
model [0. Unfortunately, for the model studied here, we can not find analytical solutions 
for these differential operators. However, we can perform an approximate analysis, and, in 
particular we can find explicitly the negatives and zero modes for the differential operator 
for the bounce wall field configuration. By making use of the spherical symmetry of the 
bounce solution, we can express the eigenvalue equation: 



in the form 



[-D E + f/"(0 wall (r - R))]Mr, 0, <p) = K*V)Mr, 0, <p) 



d 2 2d 1(1 + 1) 2 A , 2 / „ N cr . 

•3-2 - "3- + ^-T 1 + m + o^wall(^ ~R) + 7T0tall (X ~ R) 

dr z r dr r z 2 4! 



^(r) = El &ll (i)A{r) 



(37) 



where I — 0, 1, 2, Making ipi(r) = Xi( r )/ r an d z = r — R, we obtain 



d 2 1(1 + 1) 2 A , o / \ cr , 4 . . 
d? + ITTW + m+ 2^ {Z) + 4!^ ll(2) 



Xn,l{ Z ) = E Ln( n , l )Xn,l{z) ■ (38) 



Since AR <C R, we can take 1(1 + l)/(z + R) 2 « 1(1 + 1)/R 2 and then 

d 2 A 



(/ . 2 + T^wall(^) + ^tallW 



Xn(*0 = VnXn(z) 



(39) 



where r? n is obtained from 
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2 2 Kl + l) 
Vn + m 2 + -^r 1 



(40) 



We know that the — O e + U"(4>b(r)) operator has three zero eigenvalues coming from 
the bounce translational invariance. Then, for I = 1 and to the lowest value of i] n (which 
can be chosen as 771), we will have -E^ all (l, 1) = 0, with multiplicity three, as expected, 
and rjf = —m 2 — 2/R 2 . The lowest eigenvalue E 2 all (the negative eigenvalue) will be 
-2/R 2 , with multiplicity one, just what one would expect for the metastable 



state, the existence of only one negative eigenvalue |L6| . To evaluate the other eigenvalues, 
we make the following change of variable w = \f§4>o z and use (|T0f) in ([E]). We then get 



24 



dw 2 1 + e^ u 



+ 



30 



(1 + 



,V8v 



Xn{w) = V 2 n Xn{w) 



(41) 



where u 2 = r/ 2 6!/(<t0q). We can then express the eigenvalues of —D E + {7"(0 wa u(r — Rj) 
as 



-^wall 



(n, I) = 0o^[ v l + m2 + 



1(1 



R 2 



(42) 



We were not able to find any analytical solution for Eq. (0), which it is even harder 
to solve, due to the boundary conditions. We are currently working on the numerical 
solution for the eigenvalues, whose results will be reported elsewhere. 



4 The Nucleation Rate at Finite Temperature 

At finite temperature, the bounce solution is a static solution of the field equation of 
motion, 55^/5^ |^=0 B (p) = 0, where 4>b(p) is given as in fljjD , with bubble radius p and 
thickness Ap p. wall = wa ii(p) is still expressed as in (|T3D. To calculate the determi- 
nantal prefactor in we consider the eigenvalue equations for the differential operators: 

[-O e + U"(<j> B )\,l) B (i) = ti%if> B (i) (43) 

and 
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-□l 



+ U"{<j> f )]^ f {i)=^ f {i) 



(44) 



where, in momentum space, /i 2 = u 2 + E 2 , where 



27m 



are the Matsubara frequencies 



n 



0, ±1, ±2, . . .). Using ([43"D and (43) in (ph, we obtain for K the expression: 



det[-n E + U"(<p B )}p 
det [-n^ + C/"^)^ 



-1/2 



exp 



exp 



exp 



In 



det [-U E + U"{(t> B )\ 
det [-D E + U»(<f> f )h 



In 



In 





-oo Ili 






oo rij 






oo ( W n 





(45) 



where we have separated the negative and zero eigenvalues in the numerator of Eq. ([45|), 
with the prime meaning that the single negative eigenvalue, E^_, and the two zero eigen- 
values, E 2 (related to the now two-dimensional space), were excluded from the product. 
The term for n = in (uj 2 + Eq) can be handled using collective coordinates method, as 



in [10], |T5||, resulting in the factor V2 



AE 



2vrT 



where V2 is the "volume" of the two-space. 



Separating the n = modes from both the numerator and denominator of (flop and using 
the identity 



we get 



K/3 = V 2 





n 




n=l 


-AE- 




[2nT_ 


exp | 



n jT 00 ( 1 a 2 \ _ sinh(Tra) 
n=i \ n 2 ) vra 



(46) 



exp < - ln(£'2)2 - In 



sin(f|^|) 



2 \- LJ - 



+ a + b' 



(47) 



where 



-3+E-E' E'-E m 



n=l 



(4? 
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and 



E f (j)+\n(l-e-P E f 



E B (i) +\n(l-e-P EB{l) 



(49) 



Since Y,i ' has three eigenvalues less than J2j an d (E'i) 1 ^ 2 = i\E_\, we obtain 



K R = -iVo 



AE 
2^f 



2/3 2 sin (£\E- 



(50) 



Using the above results in we then obtain for the nucleation rate, per unit volume, 
the expression: 



Vi 



with 



Q 



2QT 3 exp 



/ AE \ 



AF(T) 
T 



(51) 



\E-\ 
2T 



\2nTj vrsin(^) 



(52) 



and 



AF(T) = AE-Y, 



I^ (i ) + iln(l-e-W) 



is B (i) + ilB(i-e-^« 



(53) 

In the thin-wall approximation (in an analogous way as was done in Sec. Ill) we can 
write 



AF(T) ~ -Trp 2 AC/^ (T) + 2tt(S + S p )p 
where AU e s(T) is given by 



(54) 



AU eS (T) = e + J 



d 2 p 



1 _ e -l3^/p 2 +U"(ct> f ) 



-T 



d 2 p 

(2tt) 2 



In 



1-e 



-Py/^+U"^) 



(55) 
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and Sp is the 1-loop finite temperature correction to the bubble surface energy density, 
given by 



T 

2np 



E ' 



) + In (l - 



E 



E f (j)+\n(l-e 



(56) 

where E^ all (i) are the eigenvalues of —V 2 + U"((f> wa n(p — p)). In (0), the first integral 
is divergent but it can be handled just in the same way as in the previous section, by 
the introduction of the appropriated counterterms of renormalization. The remaining 
integrals in fl55"D are all finite and they can be reduced to integrals of the type 



I(t) 



dx 



x 



ln(l 



— e 



(57) 



and it is evaluated in Appendix B. The result is: 



. . t t t 1+2? -1 " +1 C (2n) , , 

Ut) = 1(0) + - + n t__y \ > ) y > t 2 

yj yj 6 4 2 2 ^ n n + 1 2vr ) 2n 



n+l 



(5f 



where £(n) is the Riemman Zeta function. Using (|58|) in (|55|) , we obtain for (54) the 
following expression 



AF(T) = -up 2 [e + c(0/) - c(&)] + 2vrp [«5 + Sp\ , 
where (using t = f3U" \<p) 1/2 ) 



(59) 



(60) 



The critical radius for bubble nucleation, p c , is obtained by minimizing (|59|), 
5AF(T)/5p\- p=Pc = 0. 

For the eigenvalues i£ 2 all (i) of —V 2 + U"((j) wa ii(p — p)), using the now cylindrical sym- 
metry of 0waii(p) ; the eigenvalue equation 

[-V 2 + t/"(0 wall (p - p))]9i(p, <p) = <P) 



can be written as 
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d 



Id 



A 



ap z p p p z 2 4! 



where s = 0, ±1, ±2, . . .. Taking ift(p) = x(p)/p^ 2 , we obtain 



(61) 



0T 4S 2 - 1 2 A , o , , CT , 4 . _, 

rf^ + +m + 2 ^ll(P - P) + ^0wall(P " P) 



Xn,s{p) = E Lll(n, s)Xu,s(p) ■ 



As in the previous section, we make z = (p — p) and because Ap <C p, then 



<7 



'^2 + 2^ wa11 ^'' + 4j^waUV 



Xnfc) = P^Xn(p) 



(62) 



(63) 



where p 2 is now obtained from 



^Lu(n,s) =p 2 +m 2 + 



4s 2 -1 



Ap 



7-a 



(64) 



Analogously to the zero temperature case, the differential operator — V 2 + U" {(^b) has now 
two zero eigenvalues (related to the translational modes in the two-dimensional space). 
The multiplicity of fl6"4|) is two for all s ^ and for s = the multiplicity is one. 
Then for s — 1 and lower rj n (we choose pi) we will have -E^ all (l, 1) = and then 
rjl = — [m 2 + 3/(4p 2 )]. For the negative eigenvalue we obtain E 2 = E 2 a]1 (l,0) = —1/p 2 . 



As in the previous section, by taking w = Jcr/6\ <pl z, we obtain for the eigenvalues 



El a it(n, s) = 0^ v 2 n + m 2 + 4 
where v n is the same as in the previous section. 



Ap 



(65) 



5 Conclusions 

In this paper we have studied the evaluation of the vacuum decay rates, at both zero 
temperature and at finite temperature, for the cr0 6 model in D = 3, when the parameters 
of the model satisfies the conditions given in ([!]). Our main results were the determi- 



nation of the expression for the bounce solution, Eqs. ([12]) and flT9|), and also, despite 
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of the difficulties for finding the solutions for the bounce's wall eigenvalue problem, by 
taking consistent considerations for the field equations we have at hand (like the thin-wall 
approximation), we were able to perform a detailed analysis of the bounce negative and 
zero eigenvalues. We have also given a set of eigenvalue equations, which can be use- 
ful in a more detailed analysis of this problem using, e.g. numerical methods. In jy] it 
was analyzed the phase structure for the a<f> 6 model in D = 2 in the lattice and it was 
also studied the possible production of topological and nontopological excitations in the 
model. By remembering that at high temperatures our model in D = 3 resembles the 
D = 2 model at zero temperature, it will be interesting to apply the method and results 
we have obtained here to the problem studied in [p]] for the regime of a first order phase 
transition. 
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Appendix 

We evaluate here the integral fl57|), 

I(t) 

Using that 



roc 

/ dx 


xln(l-e Vx ' 2+t ' 2 ) 


Jo 





(A.66) 



dl(t) 
d(t 2 ) 



dx 







x 



d(t 2 



In ( 1 - e x 



dx 







x 



d(x 2 ) 



In ( 1 - 



(A.67) 



and integrating by parts, we obtain that 



d(t 2 ) 2 V 



-llnfl- 

9 V 



(A.68) 
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We can then write I(t) in the form: 



I(t) = 1(0) - J* dt [t In (l - e - ')] = 1(0) - J* dt 



t \ t 2 
ilnsinh ( -J + tin 2 - — 



(A.69) 



If one uses ( ff5| ) and expanding the logarithm in the above equation, we are able to perform 
the integral and we obtain the result shown in Sec. IV, Eq. (|58|). 
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Figure Captions 

Figure l:The potential U(<f>) for parameters m,\,a satisfying Eq. (1.1). 
Figure 2:The inverted potential —U{4>). 
Figure 3:The bubble field configuration wa n(r). 
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